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PThis  paper  presents  the  ''esults  of  recent  tropical  cyclone 
research  at  Colorado  State  University  (CSU).  Particular 
attention  is  paid  to  new  findings  which  impact  storm  fore- 
casting and  modeling  efforts.  Observational  studies  using 
large  amounts  of  compos ited  raw insonde,  satellite  and  aircraft 
flight  data  have  been  performed  to  analyze  tropical  cyclone 
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motion.  The  genesis  of  storms  from  earliest  stage  cloud 
clusters  was  also  analyzed.  Numerical  modeling  studies  on 
cloud  cluster  dynamics  and  the  intensification  of  pre-storm 
clusters  are  underway.  Future  research  in  all  of  the  above 
areas  is  planned  to  exploit  fully  the  very  large  tropical 
cyclone  data  set  which  has  been  assembled  at  CSU. /v 
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1 . INTRODUCTION 

The  authors  have  undertaken  a comprehensive  study  of  the  origins 
and  characteristics  of  tropical  cyclones.  Many  of  the  current  results 
are  believed  to  be  relevant  to  forecasting  and  numerical  simulations  of 
these  storms,  and  this  report  is  issued  to  disseminate  such  information. 
The  research  project  includes  composite  studies  of  large  amounts  of 
rawinsonde  and  aircraft  flight  data  and  analysis  of  digitized  DMSP 
satellite  data. 

Chapters  3-6  of  this  report  discuss  results  applicable  to  the 
forecasting  of  tropical  cyclone  movement,  intensity  changes,  genesis 
and  intensity.  Chapter  7 presents  new  information  concerning  cyclone 
structure  and  dynamics  which  are  of  particular  interest  to  numerical 
modeling  studies.  There  are  more  detailed  project  reports  (P.I.  - 
Wm.  M.  Gray)  on  each  of  the  above  topics  available  or  forthcoming  from 
Colorado  State  University  (Shea,  1972,  Shea  and  Gray,  1973  and  Gray  and 
Shea,  1973;  Ruprecht  and  Gray,  1976;  Gray,  1975a,  1975b;  George,  1975; 
Frank,  1976;  and  Zehr,  1976).  Three  other  reports  are  planned  within 
the  next  year. 
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2.  DATA  AND  ANALYSIS  TECHNIQUES 


Data  Set.  The  Information  In  this  paper  is  taken  from  various 
studies  of  tropical  cyclones  using  composited  riiwinsonde  data,  climato- 
logical information,  NOAA  flight  data  and  DMSP  satellite  photographs. 

The  data  sources  used  are  listed  below: 

1)  Ten  years  (1961-1970)  of  N.W.  Pacific  rawinsonde  data 
(*^18,000  soundings)  from  30  stations  as  shown  in  Fig.  1. 

This  data  sample  is  currently  being  expanded  to  20  years. 

2)  Twenty  years  (1956-1975)  of  N.  Atlantic  rawinsonde  data 
from  the  stations  shown  in  Fig.  2. 

3)  All  available  reduced  NOAA  Research  Flight  Facility  (RFF)  Atlantic 
inner-hurricane  flight  data  for  the  period  of  1957-1967.  This 
storm  flight  data  has  been  described  by  the  author  and  his 
graduate  student  in  published  papers  (Shea  and  Gray,  1973; 

Gray  and  Shea,  1973).  This  aircraft  data  sample  includes 
over  500  individual  radial  leg  missions.  The  data  Itself 
is  listed  in  a report  by  Gray  and  Shea  (1976).  This  is 
believed  to  be  the  best  available  information  on  inner 
tropical  storm  winds,  radius  of  maximum  winds,  inner  radar 
eye  radius,  central  pressure,  etc.  It  is  hoped  that  uses 
can  be  found  for  this  in  the  Pacific  where  some  degree  of 
meshing  of  the  aircraft  data  with  the  rawinsonde  data  may 
be  possible. 

A)  Five  years  (1971-1975)  of  direct  read-out  from  Guam  Defense  Meteor- 
ological Satellite  Program  (DMSP)  satellite  photographs  of  tropical 
cyclones  and  cloud  clusters  in  the  N.W.  Pacific.  This  1/3  nau- 
tical mile  resolution  data  has  been  digitized  and  composited  to 
perform  quantitative  analyses  of  the  convection  associated  with 
West  Pacific  tropical  weather  systems. 

5)  All  of  the  Joint  T5rphoon  Warning  Center  (JTWC) , Guam  typhoon  sum- 
maries for  the  30-year  period  of  1946-1975. 

6)  Seasonal  sea-surface  temperature  and  thermocllne  data  in  5° 
marsden  squares  for  the  whole  Pacific  as  recently  furnished 
by  the  Navy  Oceanographic  Office. 

7)  U.S.  Navy  data  tape  documentation  of  Pacific  typhoons  and 
tropical  storms  for  the  period  1946-1975. 

8)  Digitized  daily  satellite  data  for  the  period  of  October  1966 
through  1970. 

9)  The  new  global  tropical  storm  data  tape  which  has  recently 
been  issued  from  Asheville  by  Harold  Crutcher. 


Fig.  1.  Northwest  Pacific  rawinsonde  data  network. 

10)  West  Indies  tropical  cyclone  track  data  as  supplied  by 
Charles  A.  Neumann  of  NOAA  Miami. 

Rawinsonde  Compositing  Philosophy  and  Technique . 


^ Tropical  cyclones  and  cloud  clusters  spend  most  of  their  lifetimes 

over  the  warm  tropical  oceans.  Traditional  data  sources  are  very  sparse 
In  such  regions,  and  dally  tropical  weather  analyses  are  notoriously 
unreliable.  The  severe  winds  found  in  tropical  cyclones  further  reduce 
the  availability  of  such  data.  It  is  not  possible  to  obtain  enough 

I rawinsonde  data  or  surface  observations  around  any  Individual  storm  or 

cluster  at  one  time  period  to  permit  quantitative  analysis  of  structure, 

I 

* dynamics  or  energetics.  High  resolution  satellite  data  are  capable  of 

achieving  the  necessary  visual  data  density,  but  the  picture  data  are 
difficult  to  quantify,  and  the  number  of  parameters  which  can  be 
measured  accurately  Is  Inadequate  at  this  time.  In  addition,  vertical 
resolution  is  often  limited  by  a dense  cirrus  shield  which  may  cover 

[ 

most  of  the  active  convective  area. 

Aircraft  data  have  provided  the  best  Information  concerning  the 
activities  in  the  intense  central  core  regions  of  tropical  cyclones. 
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Fig.  2.  North  Atlantic  Rawinsonde  Data  Network. 

There  are  a number  of  case  studies  of  Individual  storms  based  on  Northwest 
Atlantic  hurricane  flight  data  (Rlehl  and  Malkus,  1961;  Miller,  1962; 

Gray,  1962,  1967;  LaSeur  and  Hawkins,  1963;  Sheets,  1967a,  1967b,  1968; 
Hawkins  and  Rubsam,  1968;  Hawkins  and  Imbembo,  1976)  and  also  statis- 
tical treatments  of  the  flight  data  (Shea  and  Gray,  1973,  Gray  and 
Shea,  1973).  However,  logistical  considerations  have  limited  the  ability 


of  aircraft  to  provide  Information  concerning  the  outer  convective 
regions  of  the  storm  and  its  broader  scale  environment.  Aircraft  data 
also  have  been  limited  to  a few  flight  levels  per  storm  time  period  due 
to  the  usually  low  number  of  available  aircraft,  maximum  aircraft  ceilings 
of  200  mb  or  less,  and  dangerous  low  level  flight  conditions.  These 
latter  two  restrictions  are  particularly  Important  since  the  present 
study  shows  that  the  maximum  inflow  and  outflow  from  Northwest  Pacific 
typhoons  usually  occur  near  cloud  base  and  150  mb  respectively.  Few 
aircraft  studies  of  tropical  clovid  clusters  were  accomplished  prior  to 
GATE. 

None  of  the  above  data  sources  can  produce  an  accurate  vertical 
profile  of  the  radial  wind  pattern  around  a system.  Without  such  a 
profile  it  is  Impossible  to  compute  meaningful  budgets  of  energy,  water 
vapor,  momentum,  vorticity,  etc.  In  addition,  vertical  profiles  of  the 
other  dynamic  and  thermodynamic  variables  cannot  be  determined  fully 
over  the  mesoscale  area.  It  is  necessary  to  composite  very  large 
amounts  of  data  from  many  similar  weather  systems  at  many  time  periods 
to  obtain  meaningful  quantitative  measurements. 

Although  the  extreme  variabilities  and  individual  asymmetries  of 
tropical  cyclones  and  cloud  clusters  are  well  known,  the  natures  of  the 
basic  dynamic  and  energetic  processes  which  govern  these  systems  must 
be  largely  invariant.  Compositing  allows  quantiatlve  analyses  of  these 
features.  Any  compositing  system  smoothes  out  many  of  the  individual 
characteristics  of  single  systems,  but  a great  deal  of  Information 
concerning  asymmetrical  or  "eddy"  qualities  can  be  deduced  by  the  use 
of  proper  data  handling  techniques.  The  methods  used  in  this  study 
are  discussed  below.  Storms  arc  so  variable  at  individual  time  periods 
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that  to  form  a general  physical  picture  of  their  nature.  It  Is  often 
desirable  to  smooth  out  Individual  day  variability. 

Compositing  was  performed  on  a 15°  latitude  radius  cylindrical 
grid  extending  from  sea  level  to  50  mb.  The  system  circulation  center 
was  located  at  each  time  period  using  JTWC  reports  and/or  satellite 
photographs,  and  the  grid  was  positioned  with  the  system  at  grid 
center  of  the  lowest  level.  Whenever  available  rawlnsonde  soundings  fell 
on  the  grid  at  a given  time  period  for  a given  storm,  each  sounding  was 
located  relative  to  the  storm  center  in  cylindrical  coordinates.  Figure 
3 shows  the  grid  and  the  number  of  soundings  per  grid  space  for  a typi- 
cal stratification.  All  of  the  parameters  to  be  composited,  whether 
directly  measured  or  computed  from  the  directly  measured  parameters, 
were  determined  at  the  observation  station  locations  at  19  vertical 
pressure  levels. 

The  geographical  alignment  of  the  grid  varied  with  the  coordinate 
system  used.  After  all  parameters  were  either  measured  or  computed  for 
each  sounding,  the  value  of  each  parameter  was  assigned  to  a point  at 
the  center  of  the  grid  box  in  which  the  sounding  fell.  All  soundings 
which  fell  in  that  grid  space  for  the  particular  group  of  storms  and 
time  periods  being  analyzed  were  composited. 

The  data  set  was  sufficiently  large  to  allow  compositing  of 
various  subsets.  Data  could  be  grouped  according  to  any  characteris- 
tics observed  in  individual  systems  such  as  location,  season,  intensity, 
motion,  or  intensity  tendency.  By  comparing  the  composites  of  dif- 
ferent types  of  systems  it  was  possible  to  quantitatively  analyze  the 
persistent  differences  between  the  groups.  It  was  also  possible  to 
remove  obviously  atypical  systems  or  time  periods  from  a data  group  to 
improve  the  quality  of  the  data  set. 
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The  relative  positions  of  the  system  and  the  balloon  changed  due 
to  their  respective  motions  during  the  balloon’s  ascent  time.  These 
motions  were  estimated  from  the  data,  and  the  positions  of  each  were  cor- 
rected at  each  pressure  level.  In  this  study  horizontal  eddy  fluxes  were 
estimated  by  compositing  individual  fluxes  of  quantities  and  comparing 
them  to  mean  fluxes.  The  radial  winds  were  initially  composited  and 
mass  balanced  from  the  surface  to  100  mb  by  adding  a small  constant 
correction  factor  (AV^)  to  each  individual  radial  wind  value  in  a given 
radial  band.  Changes  in  mass  of  the  volume  within  the  radial  band  were 
neglected.  For  each  sounding  the  product  of  the  corrected  and  the 
quantity  being  analyzed  was  computed  at  each  level.  Such  products  for 
all  of  the  soundings  in  an  octant  were  then  composited  as  before, 
giving  a mean  transport  value  for  each  octant  at  each  level,  V^Q,  where 
the  bar  denotes  time  and  space  averaging  of  the  V^‘Q  products.  By 
subtracting  the  product  of  the  mean  and  the  mean  quantity  (Q)  one 
could  achieve  a good  estimate  of  horizontal  eddy  transport: 

V 'Q’:^V  -Q  - iT- Q 

Characteristics  of  DMSP  Satellite  Data  and  Reduction  Procedures. 

The  data  used  were  positive  transparency  Defense  Meteorological  Satellite 
Program  (DMSP)  satellite  photos  obtained  from  the  USAF  DMSP  site  located 
on  the  Island  of  Guam  in  the  western  North  Pacific  through  the  Satellite 
Archive  Center  of  the  Space  Science  Center  of  the  U.  of  Wisconsin.  This 
data  covered  the  period  from  1971  to  1975  and  included  as  much  of  the 
life  cycle  of  each  tropical  cyclone  as  could  be  obtained  by  the  Guam 
site.  This  Included  approximately  four  days  of  data  prior  to  the  time  the 
tropical  cyclone  was  classified  as  a tropical  depression.  In  addition  150 
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photographs  of  non-developing  cloud  clusters  were  obtained  to  compare 
to  pre-storm  clusters. 

For  a detailed  satellite  study  of  penetrative  convection,  it  is 
absolutely  necessary  that  one  can  count  the  smallest  possible  diameter 
cells  which  have  penetrated  the  cirrus  canopy.  The  DMSP  VHR  imagery 
from  Guam  represents  a unique  data  set  from  the  standpoint  of  overall 
quality,  resolution  (1/3  n.ml.),  and  the  number  of  tropical  cyclones 
that  have  been  observed  over  a relatively  long  and  continuous  period 
by  the  same  type  of  satellite  sensor. 

The  Colorado  State  University,  Department  of  Atmo  .)herlc 

3 

Science's  new  Optical  Data  Digitizer  Display  (abbreviated  OD  ) is  a 
specially  designed  research  oriented  system  which  provides  the  capa- 
bility for  digitizing  transparencies.  As  many  as  16  separate  brightness 

levels  can  be  distinguished  with  the  system.  The  primary  purpose  of  the 
3 

OD  is  to  allow  the  researcher  to  selectively  prepare  portions  of  data 
stored  on  film  for  further,  more  extensive  processing  in  the  large 
university  computer. 

Deep  penetrative  convection  cells  within  a tropical  disturbance 
or  tropical  cyclone  display  an  extremely  bright  albedo,  probably 
greater  than  90  percent.  Since  this  albedo  far  exceeds  that  of  lower 
level  clouds  and  the  sea  surface,  it  is  possible  to  isolate  these  ills 

3 

using  the  OD  . The  specified  density  range  mentioned  above  was  deter- 

3 

mined  by  the  OD  as  part  of  the  brightness  set-up  for  each  new  trans- 
parency. No  brightness  values  examined  on  any  data  thus  far  have  been 
found  which  exceeded  those  of  the  penetrative  convective  cells. 

Only  those  grid  squares  which  contained  deep  convective  cells 
that  were  part  of  the  storm  system  were  dlcltized. 


Cells  contained 
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within  a erid  square  to  be  digitized  that  were  not  considered  a part 
of  the  storm  were  blocked  out. 

The  visual  DMSP  data  are  also  being  hand  analyzed  for  penetrative 
convective  elements.  The  penetrative  cloud  tops  are  outlined,  and  their 

3 

total  area  is  composited.  This  serves  as  a check  for  the  OD  since 
the  latter  system  often  has  difficulty  distinguishing  all  of  the  deep 
convection. 

Infrared  (IR)  data  with  2 n.ml.  resolution  is  analyzed  using  the 
3 

OD  to  measure  the  amount  of  high  cirrus  present  in  each  system.  The 
brightest  (coldest)  cirrus  is  stratified  into  three  levels  (temperature 
levels).  These  data  are  being  used  to  study  inter-storm  variabilities, 
diurnal  variations  in  cirrus  cover  and  lag  times  between  maximum 
penetrative  convection  and  maximum  cirrus. 

Numerical  Modeling.  We  presently  have  two  operating  6 level  PE 
numerical  models  which  have  been  developed  over  the  last  three  and 
one-half  years  by  W.  Flngerhut . One  is  3-D  and  one  axisjmmetric . The 
advectlve  terms  are  in  flux  form  and  an  eddy  diffusion  is  utilized.  The 
generalized  sigma  is  used  as  the  vertical  coorJinate.  The  tropopause, 
top  of  the  boundary  layer,  and  the  earth's  surface  are  all  coordinate 
surfaces  as  shown  in  Pig.  4.  One  can  visualize  this  sytem  as  three 
models  (stratosphere,  troposphere,  and  boundary  layer)  where  the 
boundary  conditions  are  matched. 

The  upper-most  layer  is  thought  of  as  a model  stratosphere.  It 
is  stable  and,  therefore,  resists  deformation.  This  effect  is  very 
Important  if  upper-tropospheric  and  surface  pressure  gradients  are  to 
respond  realistically  to  tropospheric  heat  sources. 


<p,  TTO  = (i)  = o 80  mb  -2 

<p,  P,  TO  = o 120  -1 

(p,  TTO  288  -.8 

(j),  TTO  456  -.6 

— Ttu,  Tiv,  tiT,  irq,  to  

(J),  uo  624  -.4 

(J),  TTO 792  -.2 

((),  P,  TTO  960  0 

(()  = O,  P,  TTO  = o 1010  1 

Fig.  4.  Vertical  structure  of  the  model  and  the  vertical  distribution 
of  variables.  Variables  at  half  levels  are  only  Illustrated 
once. 

The  boundary  layer  is  treated  as  a mixed  layer.  A frictional 
stress  is  computed  at  the  earth's  surface  and  assumed  to  decrease 
linearly  to  zero  at  the  top  of  the  layer. 

Centered  differences  are  used  both  in  space  and  time.  One  time- 
step  of  the  3-D  model  requires  about  .7  seconds  on  the  CDC  7600.  The 
2-D  model  uses  about  .02  seconds  per  time  step. 

We  have  many  experiments  planned  focussing  on  the  dynamics  of 
tropical  weather  systems,  particularly  the  cloud  cluster.  We  feel 
that  we  are  in  a very  favorable  position  to  perform  realistic  numerical 
simulations  due  to  the  availability  of  detailed  observations  obtained 


12 


by  the  weather  system  compositing  research  studies  referred  to  in 
Table  1.  We  are  planning  to  Investigate  the  topics  of  tropical  weather 
system  maintenance,  diurnal  variation,  and  tropical  storm  genesis. 

The  latter  subject  has  yet  to  be  properly  treated  in  numerical  model 
simulations.  All  previous  studies  have  dealt  with  the  intensification 
of  an  existing  closed  vortex,  with  very  large  low-level  vorticity  values, 
into  a tropical  cyclone.  Such  vortices  normally  intensify  when  located 
over  warm  tropical  waters.  The  crucial  question  is  how  these  relatively 
rare  vortices  form  initially.  More  complete  discussions  of  this  topic 
are  presented  in  Gray  (1975a)  and  Zehr  (1976). 

A crucial  step  in  the  first  area  of  investigation  was  the  diagnosis 
of  the  cumulus  forcing  in  the  cluster.  We  define  this  as  the  sum  of 
the  vertical  advectlon,  condensation,  and  evaporation  terms  in  the 
temperature  equation.  Unlike  nearly  all  existing  investigators,  we 
refer  to  the  vertical  temperature  advectlon  as  a cloud  term  since 
nearly  all  the  upward  vertical  motion  in  convective  areas  takes 
place  in  clouds  and  not  in  the  environment.  Further,  this  term  is 
very  large  as  is  the  sum  of  the  condensation  and  evaporation  terms. 
However,  the  sum  of  all  cloud  terms  is  very  small  and,  therefore, 
easier  to  deal  with.  Condensation  energy  release  is  directly  tied  to 
upward  vertical  motion.  It  typically  does  not  bring  about  direct 
environmental  sensible  temperature  gains  but  goes  primarily  into 
potential  energy  gains  (Gray,  1973).  The  rroposphere  experiences 
warming  due  to  condensation  only  as  a consequence  of  compensating 
dry  adiabatic  sinking  motion  (Lopez,  1973).  Much  of  this  compensating 
sinking  motion  occurs  at  locations  removed  from  the  cloud  areas.  The 
subsidence  warming  occurring  within  the  cloud  areas  is  largely  expended 


13 


CO 

o 


c 

0) 

u 

& 


o 

a> 

ON 

*3 

rH 

CO 

CO 

CO 

0 

c 

1 

oo 

bC 

4^ 

4J 

o 

rH 

c 

c 

CO 

CO 

CO 

v£> 

•H 

•o 

z 

B 

CO 

c 

O' 

TJ 

*r4 

0) 

•H 

c 

C 

U 

i-H 

u 

3 ^ 

3 /--N 

CJ 

CO 

O O 

O O 

bO 

3 

oc 

CO 

to 

CO 

u 

•H 

(0 

o 

bO 

a> 

0^ 

iH 

rs. 

3 

to 

o 

c 

O r-l 

o ^ 

w 

U-i 

v£) 

o 

o 

•H 

O 1 

O 1 

1 

z 

1-^ 

o 

T3 

o ^ 

O r-l 

CO 

< 

a> 

CO 

C 

• sD 

p.  \Q 

4^ 

< 

m 

£ 

u 

00 

3 

00  On 

O ON 

CO 

o 

ON 

3 

•H 

f-* 

O 

f-H 

CN  i-H 

Q 

z 

z 

z 

Z 

C>0 

CO 

U 

O 

•H 

m 

u 

pH 

•H 

1 

U 

v£> 

CO 

CO 

ir> 

oc 

o 

pH 

c^ 

c 

•H 

CO 

• 

•3 

Q) 

B 

c 

00 

o 

3 

0) 

z 

u 

o 

o 

z 

O 

pc: 

t> 

U 

o 

*H 

CO 

•H 

•H 

•H 

*3 

Ug 

CH 

CM 

0) 

4J 

4^ 

4J 

C 

•H 

•H 

•H 

CO 

in 

c 

c 

c 

c 

CO 

3 

U 

U 

4J 

oo 

r». 

o 

CO 

CO 

CO 

CO 

3 

CO 

01 

c 

ON 

c 

rH 

f-H 

Ph 

Pm 

Pu 

CO 

•H 

iH 

u 

o 

4-» 

4^ 

CO 

'O 

z 

1 

•H 

< 

< 

< 

• 

• 

• 

CO 

c 

sC 

U 

OO 

o 

4J 

3 

in 

(U 

• 

• 

• 

rH 

• 

3 

O 

ON 

i-H 

p^ 

Z 

z 

z 

o 

z 

Z 

Z 

a 

CO 

iH 

P. 

00 

^£> 

u 

O 

1 

3 

ON 

P 

f-H 

3 

5 

H 

3 

J 

Wi 

Z 

rH 

CO 

3 

3 

H 

4J 

3 

m 

•H 

4^ 

C 

•H 

c 

3 

•H 

On 

CO 

CJ 

o 

CO 

a> 

*3 

rH 

3 

3 

•H 

3 

3 

3 

< 

3 

C 

»-• 

4H 

C 

< 

pc: 

3 

3 

4J 

O 

3 

3 

3 

3 

O 

P 

O 

CO 

z 

CJ 

U 

4J 

4J 

z 

3 

O 

u 

u 

•H 

3 

3 

3 

3 

•3 

3 

3 

U-l 

O 

4J 

P 

P 

c 

P 

C 

Vm 

o 

u 

3 

O 

0 CO 

O 

O 

3 

4J 

4^ 

1 

M 

rH 

rH  U 

fH 

CO 

CO 

>s 

M 

4J 

U 

U hH 

U 

3 

•3 

3 

3 

>x 

>^  4J 

>x 

3 

3 

3 

3 

3 

3 

3 

C 

C 

U 

O 3 

CJ 

CJ 

•H 

•H 

•H 

4^ 

ki 

g 

C 

3 

OC 

4J 

4J 

4-t 

3 

O 

0 

g 

M 

P 

rH 

rH  U 

^H 

rH 

P 

P 

P 

rH 

4J 

U 

CJ 

3 

3 

3 3 

3 

3 

3 

3 

3 

P. 

u 

CO 

3 

4J 

U 

3 C 

U 

3 

P 

rH 

rH 

rH 

6 

3 

P 

P 

Z 

•H 

*H  UJ 

•H 

•H 

O 

4J 

4H 

4H 

o 

t-l 

3 

3 

3 

3 

00  Ox 

CL 

(X 

P. 

P 

•H 

< 

< 

< 

CJ 

,P] 

C 

C 

4J 

U 

•H  N43 

o 

O 73 

o 

O 

OO 

3 

P 

C 

3 

•H 

rH  On 

V4 

u c 

)H 

u 

3 

• 

• 

• 

CO 

M 

M 

P 

Pm  rH 

H 

H 3 

H 

H 

pc 

z 

z 

z 

0) 

c 

o 


/— N 

rH 

m 

cn 

VO 

/— s 

3 

r*. 

>N 

ON 

o\ 

ON 

m 

CJ 

CO 

rH 

rH 

rH 

3 

ON 

rH  CO 

•H 

r 

r 

r 

rH 

3 3 

*o 

>N 

3 

3 

3 hH 

3 

3 

3 

3 

r 

•H  -O 

4H 

JC 

Z 

3 

m 

P 3 

3 

CO 

O 

CO 

CO 

m 

vO 

O 4H 

Wi 

ps. 

ON 

NO 

U CO 

3 

n 

•o 

TS 

'V 

ON 

ON 

rH 

H 

4J 

P 

•H 

3 

P 

c 

p 

^H 

rH 

ON 

Z 

3 

O 

CO 

4J 

3 

3 

3 

#k 

rH 

4J  3 

4H 

3 

•H 

3 

3 

3 

P M 

3 

M 

P 

rH 

3 

>s 

>x 

«k 

00 

r 

3 3 

3 

4J 

O 

3 

a 

3 

3 

3 

c 

U 

V* 

U 3 

H-: 

CO 

z 

CJ 

B 

Z 

U 

U 

3 

3 

o 

z 

u t/i 

P 

o 

CO 

CJ 

O 

>H 

M 

3 

3 

3 3 

3 

CJ 

o 

Um 

CJ 

N3 

U P 

CO 

rH 

CM 

PO 

L 


15 


in  balancing  the  cloud  region  liquid  water  reevaporation  due  to  cloud 

detrainment.  Williams  and  Gray  (1973)  and  Ruprecht  and  Gray  (1976) 

have  demonstrated  that  the  direct  Influence  of  the  typical  4°  diameter 

2 

cluster  condensation  (equivalent  to  1000-1500  calorles/cm  per  day) 
on  cluster  sensible  temperature  change  is  nearly  zero.  Cluster 
cumulus  clouds  typically  exert  only  a very  small  sensible  heat  change 
upon  the  cluster  itself.  Most  cumulus  condensation  energy  is  exported 
to  the  surrounding  regions  to  balance  the  outer  subsidence. 

We  believe  it  is  Important  to  point  out  that  the  net  effect  of 
all  the  cloud  Induced  temperature  changes  is  much  less  than  the  clear 
atmosphere  net  radiative  cooling.  Radlatlonal  forcing  appears  to  be 
larger  than  that  of  cumulus  convection  in  the  organized  tropical  cloud 
cluster.  This  has  broad  implications  for  other  modeling  work. 

Our  work  on  the  modeling  of  the  diurnal  variation  of  the  cluster 
is  well  under  way.  We  have  reproduced  the  large  diurnal  variations  in 
rainfall  rates  and  divergence  profiles  which  were  observed  by  Jacobson 
and  Gray  (1976).  This  shows  a significant  diurnal  variation.  We  have 
also  noted  two  ways  in  which  the  dlurnally  varying  radiation  affects 
the  cluster.  The  obvious  one  is  through  the  time  varying  radiative 
cooling  of  the  surrounding  clear  environment.  However,  the  energy 
gain  during  the  day  (loss  at  night)  from  cloud  tops  in  the  cluster 
is  equally  important.  Both  Influence  the  hour  of  maximum  rainfall 
and  the  character  of  the  cluster  inflow  profile. 

Summary  of  Recent  and  Current  Research  Projects.  Table  1 shows 
the  tropical  cyclone  research  projects  which  are  either  now  in  progress 
or  have  been  completed  since  1972  at  CSU.  The  more  Important  findings 
of  these  studies  are  summarized  in  the  following  chapters.  Many  of  the 


results  presented  in  this  report  have  applicability  to  forecasting  and 
numerical  modeling. 


r 
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3.  TROPICAL  CYCLONE  MOTION 

Steering  Flow.  The  N.W.  Pacific  rawlnsonde  compositing  study  by 
George  (1975)  revealed  strong  correlations  between  the  direction  and 
speed  of  storm  motion  and  the  area  weighted  average  flow  from  1-7*^  radius. 
It  was  found  that  tropical  cyclones  move  about  16*^  to  the  left  of  the 
average  1-7°  flow  at  500  mb.  Their  speeds  of  movement  average  1.16 
times  the  speed  of  the  1-7°  wind  at  700  mb.  Other  radial  bands  and 
vertical  levels  did  not  give  such  consistent  results. 

The  above  values  are  not  merely  the  result  of  averaging  many  storms 
with  greatly  different  storm  motion-surrounding  flow  relationships.  The 
data  were  divided  into  13  stratifications  according  to  latitude,  speed 
and  direction  of  movement,  intensity  and  intensity  change.  Twelve  of 
these  stratifications  are  Independent  subsets  of  the  typhoon  data  set 
(all  soundings  with  central  pressures  P^  £ 980  mb) , and  the  remaining 

data  set  (980  < P £ 1000  mb)  Is  completely  Independent  of  the  others. 

I 

! The  stratifications  and  their  storm  motion-surrounding  flow  relation- 

ships are  summarized  In  Table  2.  The  storm  direction  varied  only  from 
12-23°  left  of  the  averaged  500  mb  wind.  Storm  speed  was  always  from 

1 1.07  to  1.2A  times  the  700  mb  winds.  These  relationships  are  remarkably 

i 

consistent  considering  the  Independence  of  the  data  subsets  and  the  wide 
variety  of  storm  types  encompassed.  Cyclones  moving  In  various  direc- 
tions and  at  different  speeds  behave  similarly. 

Steering  flow  was  also  investigated  for  winds  averaged  through 
various  vertical  levels.  Table  3 shows  the  best  results  for  the  same 
13  stratifications  shown  In  Table  2.  The  most  consistent  correlation 
with  direction  of  storm  motion  was  found  for  the  700-500  mb  layer  winds 
while  storm  speed  fit  the  100-500  mb  averaged  winds  best.  In  both  cases 


TABLE  2 


HORIZONTAL  MHI.TIPI.E  KAPIAL  HAND  SUMMARY  FOR  1-7°  BAND 


Stratification  Mean  Stoim 

U 1 *’  f 

Masan  I -7 

Sfonu 

Mean 

Mean  l”-7° 

Ratio  Storm 

Olrectlon 

Surrounding 

Opvtnt Ion 

Storm 

Surround!  ttg 

Spred/Surround Ing 

Wln<l  IHrct’tloii 

f rop) 

Spood 

Wind  Spo«'<l 

1-7*'  Wind  Sprr<l 

at  ^00  tnl> 

Surrounding 

at  700  ah 

Wind 

. -1, 

-1 . 

ClriilU- _ (»»'')  (»"'') 


Latitude 


Latitude  >20^N 

352° 

010° 

18‘’l 

5.61 

5.06 

1.11 

Latitude  <20°N 

300° 

314° 

14°L 

5.05 

4.23 

1.19 

Speed 

Slow  Storm 

Speed  0-3  ms 

338° 

001° 

23°L 

2.43 

2.19 

1 .11 

Moderate  Storm 
Speed  3-7  aa”' 

326° 

344° 

18°L 

5.19 

4.37 

1.19 

Fast  Storm  . 

Speed  >7  as 

006° 

018° 

12°L 

10.12 

8.71 

1.16 

Direction 

Direction  A 

250°  < Dir.  < 310° 

285° 

300° 

15°L 

6.16 

4.96 

1.24 

Direction  B 

310  < Dir.  ^ 350° 

324° 

341° 

17°L 

5.32 

4.39 

1.21 

Direction  C 

350°  < Dir.  < 060° 

027° 

042° 

15°L 

7.08 

6.59 

1.07 

Intenalty 

Intenalty  1 
lOOOab  > C.P.> 

980Bb 

319° 

331° 

12°L 

4.87 

4.07 

1.20 

Intensity  2 

980Bb  > C.P.  > 

950^ 

326° 

343° 

17°L 

5,03 

4.22 

1.19 

Intensity  3 

C.P,  ^ 950rt) 

319° 

341° 

22°L 

5.17 

4.34 

1.19 

Intensity  Chanae 

Deepening  Storms 

304° 

317° 

13°L 

4.89 

4.31 

1.13 

Pilling  Storms 

360° 

017° 

I7°L 

5.69 

5.13 

i.ll 

Mvsn; 

1^.4"\. 

Mean: 

1.16 

Standard  Deviation: 

1.4° 

Standard  DovlatUm: 

o.os 
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TABLE  3 


VKRIICAI.  Min.Tiri.K  KAJUAI 

HAND  SlfMMAKY 

FOR  1 - 

7"  BAND 

St  rat  Klc  at  ion 

Menu  St4>riti 
1)1  rert  Ion 

Mean  l' 

Surroiiinl  (ng 

Wind  Direct  Ion 
(700-500  mh) 

St  t)nn 

Dev  l«t Ion 

from 

.Surroun<l  Ing 
Wind 

Mean 

Storm 

Speed 

/ 

(wm 

Mean  l°-7‘’ 

Sui  rounding 
Wind  Speed 
(1000-500  mb) 
) (ms 

Ratio  Storm 

Speed /Surrounding 
I-7”  Wind  Speed 

Latitude 

Latitude  >20^N 

352° 

007° 

15”l 

5.61 

4.55 

1.23 

Latitude  <20°N 

300° 

312° 

12°L 

5.05 

4.05 

1.25 

Speed 

Slow  Storn  ^ 
Speed  0-3  ms 

338° 

359° 

21°L 

2.43 

1.91 

1.27 

Moderate  Storm 
Speed  3-7  ms  ^ 

326° 

342° 

16°L 

5.19 

4.03 

1.28 

Fast  Storm  ^ 

Speed  >7  ms” 

006° 

017° 

11°L 

10.12 

7.77 

1.30 

Direction 

Direction  A 

250  < Dir.  _<  310° 

285° 

302° 

17°L 

6.16 

4.65 

1.32 

Direction  B 

310°  ' Dir.  350° 

324° 

337° 

13°L 

5.32 

4.25 

1.25 

Direction  C 

350®  < Dir.  ^ 060® 

027° 

040° 

n°L 

7.08 

5.67 

1.25 

Intensity 

Intensity  1 

1000^  > C.P.  > 
980mb 

319° 

333° 

14°t. 

4.87 

3.71 

1.31 

Intensity  2 

980mb  ^ C.P.  > 
950mb 

326° 

344° 

1«°1. 

5.03 

3.67 

1.29 

Intensity  3 

C.P.  ^ 950mb 

319° 

340° 

21°L 

5.17 

4.03 

1.28 

Intensity  Change 

Deepening  Storms 

304° 

318° 

14°L 

4.89 

3.82 

1.28 

Filling  Storms 

360° 

019° 

19°L 

5.69 

4.55 

1.25 

Mean 

: 15.7°L 

Mean:  1.27 

Standard  Deviation 

1 3.3° 

Standard  Deviation:  0.03 
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area  weighted  averaged  winds  from  1-7°  were  found  to  give  the  most  con- 
sistent results. 

The  storm  motion-surrounding  flow  relationships  were  essentially 
similar  when  geostrophlc  winds  computed  from  the  observed  heights  were 
used  Instead  of  observed  winds.  This  Indicates  that  height  and  wind 
fields  may  be  used  Interchangeably  for  storm  motion  forecasting  according 
to  operational  convenience. 

The  relationships  between  storm  motion  and  the  winds  in  the  octant 
to  the  right  of  the  direction  of  storm  motion  were  investigated.  It  was 
found  that  the  right  octant  winds  at  850  mb  at  all  radii  from  1°  to  9° 
were  within  about  5°  of  the  mean  storm  direction  for  all  stratifications. 
There  were  no  good  correlations  between  right  octant  winds  and  storm 
speed. 

Further  steering  flow  analyses  are  being  performed  for  N.  Atlantic 
tropical  cyclones  to  verify  the  results  more  thoroughly.  A comparison 
of  right  vs.  left-turning  storms  Is  being  made  to  try  and  separate  the 
factors  affecting  storm  motion  from  mean  environmental  flow.  Early 
results  Indicate  that  there  are  some  persistent  differences  in  the  200 
mb  wind  fields.  Data  prior  to  motion  change  will  be  examined  to  determine 
If  there  Is  a lag  time  between  changes  in  the  wind  field  and  changes  in 
storm  motion.  If  such  a lag  Is  found,  observations  of  surrounding  winds 
might  allow  forecasts  of  sudden  direction  changes.  In  addition,  the 
denser  data  network  In  the  West  Indies  may  allow  studies  of  the  varia- 
bilities of  Individual  storm  motions  rather  than  groups  of  storms.  This 
was  not  possible  in  the  data-sparse  N.W.  Pacific. 

Recurvature.  The  most  crucial  problem  facing  tropical  cyclone 
forecasters  Is  the  prediction  of  recurvature.  To  analyze  the  features 


21 


of  the  flow-field  associated  with  recurvature,  twenty-one  pairs  of 
tropical  cyclones  were  selected.  Each  pair  inc Hided  one  recurving  and 
one  non-recurving  storro  where  the  tracks  prior  to  recurvature  were  within 
5°  latitude  of  each  other.  The  tracks  of  the  two  types  of  stoms  are 
shown  in  Figs.  5 and  6.  A separation  point  (S)  was  identified  where 
the  recurving  storm  began  to  deviate  significantly  from  the  non-recurving 
track.  Data  was  composited  for  each  type  of  storm  (recurving  and  non- 
recurving) for  various  time  intervals  up  to  60  hours  before  S. 

Differences  between  recurving  and  non-recurving  storms  were  observed 
in  both  the  lower  and  upper  troposphere,  but  upper  level  differences 
were  far  more  pronounced.  Figure  7 shows  composite  200  mb  winds  for 
both  types  of  storms  and  the  difference  between  the  two  types  60  hours 
prior  to  the  separation  point.  The  scale  of  the  cylindrical  grid  has 
been  altered  so  that  the  closest  plotted  point  is  8°  latitude  from  the 
storm  center.  The  westerlies  10-20°  north  of  the  storm  are  much  stronger 
and  exist  closer  to  the  storm  center  for  recurving  storms.  Note  that 
these  large  differences  occur  2h  days  before  recurvature.  As  the  storms 
approach  the  recurvature  point,  the  westerlies  get  stronger  and  closer 
to  the  center  (Figs.  8 and  9 ) . It  has  long  been  known  that  tropical 
cyclones  tend  to  recurve  under  the  Influence  of  westerly  troughs  to  the 
north,  but  this  is  the  first  quantitative  analysis  showing  a strong 
correlation  between  upper  level  winds  and  recurvature  well  in  advance  of 
the  first  change  of  direction.  Considerable  improvements  in  forecasting 
skill  may  result  from  careful  analysis  of  upper  level  winds.  Table  4 
summarizes  the  recurving/non-recurving  wind,  height  and  temperature  dif- 
ferences at  700  mb  and  200  mb  from  12  hours  to  60  hours  prior  to  the 
separation  point. 
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Fig.  6.  Actual  tracks  of  the  21  non-rec 
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200  mb  Dir.  & Speed 
NON-RECURVING  at  time 
S-60  hrs. 


200  mb  Dir.  & Speed 
RECURVING  at  time 
S-60  hrs. 


200  mb  Dir.  & Speed 
NON-RECURVING-RECURVING 
DIFFERENCE  at  time 
S-60  hrs. 


. 7.  Composite  200  mb  winds  for  non-recurvlnft  and  recurving  storms 
60  hours  prior  to  the  separation  point.  Also  shown  is  the 
difference  between  the  2 composites. 


200  mb  Dir.  & Speed 
NON-RECURVING  at  time 
S-12  hrs. 


200  mb  Dir.  & Speed 
RECURVING  at  time 
S-12  hrs. 


200  mb  Dir.  & Speed 
NON-RECURV  TNG-RECURV ING 
DIFFERENCE  at  time 
S-12  hrs. 
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I Current  and  future  recurvature  research  at  CSU  Is  centered  on  an 

I extension  of  the  above  analysis  to  North  Atlantic  storms  and  on  establish-  i 

ing  individual  case  variability.  It  is  hoped  that  an  operational  fore- 
casting scheme  will  result.  DMSP  data  are  also  being  analyzed  to  deter- 
mine any  relationships  between  storm  motion  and  convective  patterns, 
f Recommendat ion . Individual  case  statistics  using  200  mb  cyclone 

data  should  be  developed.  The  middle  tropospheric  winds  at  500  and  700 
mb  should  not  be  used  exclusively  to  forecast  recurvature . 

I Relationships  between  changes  in  storm  motion  and  changes  in  the 

i surrounding  wind  and  height  fields  must  be  explored  in  more  detail. 

I 

I It  is  extremely  important  to  determine  whether  there  are  observable 

changes  in  the  environment  prior  to  changes  in  storm  direction.  The 

individual  variabilities  of  storm  motion  with  respect  to  observed  en- 

I 

I 

vlronmental  features  should  be  analyzed  further. 

1 

i 

I 
f 
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4.  INTENSITY  CHANGES 

Little  skill  has  been  developed  in  the  forecasting  of  tropical 
cyclone  intensity  changes.  The  N.W.  Pacific  typhoon  compositing  study 
(Frank,  1976)  was  used  to  examine  persistent  differences  between  tropical 
cyclones  with  different  central  pressure  (P^,)  tendencies.  Four  data  sets 
were  assembled: 

1)  Deepening  Typhoons  - all  storms  with  P ^ 980  mb  and  with 
central  pressures  in  sustained  falling  trends, 

2)  Filling  Typhoons  - all  storms  with  P^  980  mb  and  with  the 
central  pressures  in  sustained  rising  trends, 

3)  Developing  Tropical  Storms  - all  storms  with  980  mb  < P^  ^ 1000  mb 
with  falling  central  pressures  and  which  soon  reached  typhoon 
intensity  (P^  < 980  mb) , and 

4)  Non-Developing  Tropical  Storms  - storms  with  980  mb  < P^  _<  1000  mb 
which  never  reached  typhoon  intensity. 

The  results  are  as  follows: 

Potential  Buoyancy.  The  value  of  0^,  at  the  surface  minus  the 
saturated  0 (0^c)  value  at  600  mb  gives  a measure  of  the  atmosphere's 

potential  to  sustain  deep  cumulus  buoyant  moist  convection.  Values  of 
®^sfc  ~ ®ES600  plotted  for  deepening  and  filling  typhoons  (Cases 

1 and  2)  in  Figs. 10  and  11.  Note  that  the  potential  buoyancies  to  the 
northwest  of  the  storm  centers  differ  considerably.  (The  mean  direction 
of  storm  motion  is  NNW) . Filling  storms  tend  to  have  large  areas  of 
strongly  negative  buoyant  air  in  this  region  while  deepening  storms  are 
closer  to  neutral  buoyancy  in  these  areas.  The  difference  seems  to 
arise  from  cooler,  drier  surface  air  in  the  filling  storms.  Similar 
results  are  found  for  developing  storms  (neutrally  buoyant)  and  non- 
developing  storms  (negatively  buoyant).  Since  substantial  low  level 
inflow  relative  to  the  moving  storm  center  comes  from  the  northwest,  it 
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is  clear  that  the  filling  and  non-developing  storms  are  advecting  In  low 
level  air  which  is  relatively  unfavorable  for  moist  convection.  This  may 
be  an  important  influence  on  the  weakening  or  non-developing  storms.  In- 
tensifying storms  also  show  somewhat  more  favorable  buoyancy  conditions 
to  the  south  than  do  weakening  storms . 

The  above  results  are  essentially  unchanged  when  the  data  is  stra- 
tified by  latitude  (North  and  South  of  20°N) . However,  there  are  season- 
al and  latitudinal  biases  in  the  data.  The  percentages  of  storms  occur- 
ring in  the  winter  and  spring  are  significantly  larger  for  filling  and 
non-developing  storms  than  for  deepening  and  developing  ones.  During 
these  seasons  the  sea  surface  temperatures  to  the  northwest  of  the  mean 
storm  centers  become  quite  cold.  The  filling  and  non-developing  storms 
also  tend  to  be  farther  north  than  their  Intensifying  counterparts.  It 
is  not  clear  whether  the  observed  cool,  dry  low  level  air  is  directly  re- 
sponsible for  storm  weakening.  The  fact  that  filling  and  non-developing 
storms  are  relatively  more  common  in  winter  and  at  northerly  latitudes 
suggests  a correlation,  but  other  factors  may  be  important  or  dominant. 
The  above  study  does  lend  some  support  to  the  hypothesis  that  storms  tend 
to  weaken  as  they  move  over  cold  water  as  discussed  by  Brand  (1971). 

Upper  Tropospheric  Trough.  A strong  upper  level  (300-100  mb) 
trough  has  been  observed  in  composites  of  filling  and  non-developing 
storms.  The  trough,  maximum  at  200-150  mb,  is  observed  at  large  radii 
(10-14°)  to  the  northwest  of  the  storm  center.  It  can  be  determined 
from  both  the  height  and  wind  fields.  Figures  12  and  13  show  the  150  mb 
height  and  wind  fields  for  deepening  and  filling  storms,  respectively. 

The  trough  is  cold  core.  It  is  related  to  numerous  factors  which  might 
cause  Intensity  changes  Including  the  low-level  stability  and  the 
characteristics  of  the  upper- level  outflow. 
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Outflow.  Typhoon  outflow  at  large  radii  (8-14'’)  Is  typically 
concentrated  in  channels.  Upper  level  (150  mb)  outflow  for  deepening 
typhoons  Is  maximum  to  the  southwest.  Very  little  outflow  Is  observed 
to  the  northeast.  However,  the  composite  for  filling  typhoons  shows  a 
maximum  outflow  toward  the  northeast  with  a weaker  channel  to  the  south- 
west. These  patterns  are  evident  in  Figs.  12  and  13.  When  the  data  are 
further  stratified  by  latitude  (north  and  south  of  20'’n)  the  results 
are  similar  except  that  the  southwest  outflow  chiinnel  for  filling  storms 
south  of  20”  is  stronger  than  their  northeast  channel.  Similar  but  less 
conclusive  differences  are  seen  In  the  composite  wind  fields  for  develop- 
ing vs.  non-developing  storms.  The  northeast  outflow  for  filling  storms 
is  related  to  the  composite  trough  seen  in  Fig.  13.  These  results  con- 
tradict those  of  the  recent  satellite  studies  by  Ramage  (1973)  and  Tsui 
et  al.  (1977)  which  reported  strong  cirrus  outflow  to  the  northeast  of 
deepening  storms.  This  could  result  from  a poor  correlation  between  max- 
imum outflow  at  outer  radii  with  cirrus  cover.  However,  the  Tsui  study 
showed  that  a cirrus  outflow  stream  was  correlated  with  previous  inten- 
sification, not  current  intensity  tendency,  and  Ramage  studied  only  a few 
cases.  It  is  possible  that  storms  intensify  prior  to  interaction  with  a 
westerly  trough  and  then  weaken  as  it  approaches  as  suggested  by  Riehl  (1972). 

The  physical  mechanisms  linking  outflow  characteristics  to  Intensity 
change  are  not  known.  They  are  probably  related  to  greater  angular  mo- 
mentum and  kinetic  energy  exports  in  the  filling  systems  or  to  changes 
in  the  thermodynamic  fields  such  as  extra  temperature  and  moisture 
ventilation.  Current  research  is  concentrating  on  this  subject.  Indi- 
vidual case  studies  will  be  performed  to  determine  how  well  northeast 
outflow  jets  predict  storm  weakening  on  a case  by  case  basis. 
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Convection  as  Observed  by  DMSP  Satellite  Data.  Studies  of  the 
amounts  of  deep  convection  surrounding  deepening  and  filling  storms 
using  DMSP  photographs  show  only  slight  differences.  The  deepening 
storms  tend  to  have  slightly  more  deep  convection,  but  the  differences 
are  much  too  small  to  be  useful  on  an  individual  case  basis.  Studies 
of  the  relationship  between  convective  organization  and  intensity  change 
are  planned  for  the  immediate  future. 

Summary . It  is  extremely  difficult  to  predict  storm  intensity 
changes.  The  differences  between  weakening  and  intensifying  storms  are 
small  and  are  usually  obscured  by  large  inter-storm  variabilities.  How- 
ever, a few  differing  characteristics  are  noted.  Tropical  cyclones  tend 
to  weaken  as  they  approach  low  level  cold  dry  air  - usually  associated 
with  cold  water.  In  the  mean  an  upper  level  trough  to  the  northwest  and 
an  outflow  jet  to  the  northeast  are  indicative  of  typhoon  weakening  or 
of  non-intensification  of  a tropical  storm  (980  < P^  1000  mb).  No 
obvious  correlation  between  the  amounts  of  deep  convection  and  intensity 
tendency  are  observed. 

Future  Research.  The  primary  factor  limiting  composite  studies  of 
tropical  cyclone  intensity  changes  is  an  insufficient  number  of  storm 
cases  in  each  group.  The  Pacific  data  set  is  being  expanded  to  20  years 
doubling  the  size  of  that  sample,  and  a 20  year  data  set  is  also  being 
assembled  for  the  West  Indies.  It  is  anticipated  that  the  combined  data 
from  these  two  data  sets  will  allow  more  quantitative  analyses  of  the 
processes  causing  storm  intensity  changes.  In  addition,  the  larger  data 
set  will  allow  stratification  of  storms  into  more  specific  groups. 

Storms  undergoing  intensity  change  may  be  further  stratified  by  season, 
latitude,  rate  of  central  pressure  change,  current  intensity,  etc.  This 
will  be  a prime  topic  of  future  research. 
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5.  TROPICAL  CYCLONE  GENESIS 

The  climatological  conditions  favorable  to  the  formation  of  tropical 
cyclones  have  been  summarized  by  Gray  (1968,  1975a),  but  the  physical 
processes  which  occur  are  not  well  documented.  It  is  known  that  they 
develop  from  rather  ordinary  looking  cloud  clusters  or  disturbances. 
However,  no  clear  differences  between  common  ordinary  clusters  and  rela- 
tively rare  pre-typhoon  clusters  have  been  observed. 

Zehr  (1976)  compared  10  years  (1961-1970)  of  N.W.  Pacific  rawinsonde 
data  around  clusters  which  developed  into  typhoons  with  data  computed 
from  two  years  of  clusters  which  did  not  intensify  into  cyclones.  He 
examined  ITCZ  clusters  only  (latitude  5-18°  N) . Trade  genesis  at  higher 
latitudes  was  not  treated.  Overall  the  clusters  were  very  similar,  but 
some  significant  differences  were  found. 

In  this  study  stage  0 refers  to  all  ITGZ  clusters  which  did  not 
develop  into  tropical  cyclones.  Stage  00  is  restricted  to  the  months 
of  June-September  and  to  longitudes  from  125°  E to  165°  E.  Stage  2 
refers  to  pre-typhoon  clusters  at  all  time  periods  up  to  1 day  prior 
to  the  first  aircraft  reconnaissance  observation  (usually  about  30 
knots  maximum  wind).  The  most  relevant  comparison  of  developing 
to  non-developing  clusters  is  between  stages  00  and  2 since  the  stage  00 
set  is  more  representative  of  tropical  cyclone  season  and  locations  than 
the  stage  0 set. 

The  same  subject  is  being  investigated  using  DMSP  satellite  photo- 
graphs to  analyze  convective  differences  (Erickson,  1977).  The  more 
important  results  of  both  studies  are  discussed  below. 

Divergence  and  Vertical  Motion.  Developing  clusters  (stage  2)  tend 
to  have  more  convergence/divergence  at  4°  radius  and  a greater  maximum 
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r upward  vertical  motion  (w)  than  non-developing  clusters  (stage  0 and  00), 

(rigs.  14  and  15).  The  difference  In  w Is  about  25%.  This  Implies 

max 

[ more  convection  In  the  developing  clusters  and  Is  In  good  agreement  with 

i the  DMSP  study  showing  more  deep  convection  from  0-4°  for  such  clusters. 

Most  of  the  increased  convection  observed  in  developing  clusters  occurs 
j between  2°  and  4°  radius.  There  Is  little  difference  between  the  amounts 

I 

! of  mean  convection  in  the  center  region  (0-2°) . Unfortunately,  Individual 

I clusters  are  so  variable  that  storm  genesis  is  probably  not  reliably 

j predictable  from  observations  of  amounts  of  cluster  convection, 

j Low  Level  Vortlclty.  It  Is  well  known  that  tropical  cyclones  tend 

to  form  in  regions  with  strong  low  level  vortlclty  such  as  the  ITCZ 
(Sadler,  1967;  Gray,  1975a).  The  900  mb  relative  vortlclty  fields  of 
the  mean  pre-typhoon  cluster  and  of  the  non-developing  cluster  data  sets 
are  shown  in  Fig.  16.  Note  that  the  area  with  values  greater  than  2 x 
10  ^ sec  ^ is  significantly  larger  for  the  pre-typhoon  cluster.  In 
the  immediate  vicinity  of  the  cluster,  the  relative  vortlclty  is  about 
twice  as  large  with  the  pre-typhoon  cluster  as  with  the  non-developing 
clusters.  This  is  the  most  substantial  difference  between 
ing  and  pre-typhoon  cloud  clusters.  It  may  be  possible  to 
minimum  threshold  va ' of  low  level  vortlclty  over  a 'v  4° 
for  tropical  cyclone  gen  'Sis.  Tlie  problem  Is  to  determine 
mesoscale  vortlcltles  irum  the  limited  data  available  over 
oceans.  Perhaps  satellite  det  .mined  winds  could  be  used, 
an  observed  "sharpening"  of  the  ITCZ  suggests  locally  high 
low  level  vortlclty.  The  use  of  vortlclty  observations  as 
prediction  tool  Is  currently  being  Investigated. 

, 
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Fig.  14.  Mean  d i verren . e within  the  r = 0-4*^*  area  for  the  pre-typhoon 
cluster  C"  and  the  nnn-<leveloplng  clusters  (<')  and  (00). 
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Fig.  15.  Mean  vertical  velocity  (lo)  In  the  r * 0-4^  area  for  the  pre- 
t/phoon  cluster  (2)  and  the  non-developing  clusters  (0)  and 
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Vertical  Temperature  Structure.  Developing  clusters  tend  to  be 
slightly  warm  core  in  the  lower  troposphere  even  at  the  earliest  ob- 
servation time  while  non-developers  are  cold  core  at  those  levels  (first 
noted  by  Rlehl  in  the  1940's).  Figure  17  compares  the  temperature  anomalies 
of  both  types  of  clusters.  Assuming  approximate  thermal  wind  balance,  the 
mean  tangential  wind  profiles  should  show  different  vertical  gradients  of 
V^.  This  proves  to  be  the  case  (Fig.  18).  The  dc'veloplng  clusters  have 
maximum  cyclonic  flow  in  the  lowest  levels  while  non-developing  clusters 
have  maximum  V.^,  in  the  middle  levels. 

The  small  temperature  anomalies  in  clusters  are  very  hard  to  observe 
for  Individual  cases  due  to  the  relatively  large  inter-sounding  varia- 
bilities. However,  observation  of  the  vertical  tangential  wind  shear 
holds  some  promise  as  a tropical  cyclone  genesis  forecast  aid,  especially 
since  the  warm  core  is  observed  at  the  earliest  stages  of  development. 

Anticyclonlc  Outflow.  The  200  mb  outflow  Is  more  antlcyclonic  for 
developing  than  for  non-developing  clusters  as  shown  in  Table  5.  Satellite 
determined  winds  should  be  plentiful  enough  at  this  level  to  allow  determina- 
tion of  mean  V values  at  various  radii. 

T 


TABLE  5 


Mean  Tangential  Wind  at  200  mb 
(m  sec  S 


Date  Set 

r = 2° 

r = 4 

4 = 6 

Non-developing 

-1.35 

-1.95 

-1.09 

(Stage  00) 

Non-developing 

-0.70 

-0.66 

-1.25 

(Stage  0) 

"Pre-typhoon 
(Stage  2) 


-1.18 


-1.72 


-2.90 
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Fig.  17.  Vertical  temperature  profiles:  pre-typhoon  cluster  (2)  versus 
non-developing  clusters  (0  and  00).  T^,  = mean  temperature 

of  the  cluster.  = mean  temperature  of  the  surrounding 

region.  - r = 0-1°,  : 5-7°. 


Vrlm  sec ')  r»4* 
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C (IO  *Mc')  r.0-4- 

Fig.  18.  Vertical  tangential  wind  profiles  of  the  pre-typhoon  clusters 
(2) , and  non-developing  clusters  (0  and  00) . The  mean  tan- 
gential winds  shown  are  the  mean  In  the  r * 3-5°  band,  and 
can  be  Interpreted  as  a mean  relative  vortlclty  within  r = 4°. 
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Vertical  Shear  and  Ventilation.  Strong  vertical  wind  shear  is 
unfavorable  for  storm  development  (Cray,  1973,  1968)  . Figure  19  shows 
that  the  zonal  vertical  shear  is  less  for  developing  clusters  (stage  2) 
than  for  non-developers  (stage  00).  The  weak  mean  vertical  shear  for 
stage  0 apparently  results  from  the  cancellation  of  many  strong  shears 
of  opposite  signs.  Also  shown  are  the  propagation  speeds  of  the  clusters. 

The  radial  wind  flow  relative  to  the  cluster  can  be  analyzed  in 
terms  of  its  divergent  and  non-d ivergent  ("ventllative")  components  at 
the  approximate  cluster  perimeter — r = 4^’  (Zehr,  1976).  The  divergent 
component  is  merely  the  absolute  value  of  the  mean  radial  wind  (|v^|), 
and  the  ventilation  (VEN)  is  given  hy:  VEN  = !V^|-  |v^i.  Figure  20 
shows  the  mean  divergent  and  ventllative  components  of  the  radial  wind 
for  both  types  of  clusters.  It  is  clear  that  ventilation  is  strongest  in 
the  non-developing  clusters.  Ventilation  tends  to  reduce  the  temperature 
and  moisture  anomalies  of  the  clusters. 

Summary.  Several  major  facts  concerning  differences  between  TTCZ 
pre-typhoon  cloud  clusters  and  non-developing  ITCZ  clusters  have  been 
established.  Many  of  these  reinforce  the  findings  and  conclusions  of 
previous  studies. 

1)  The  mean  vertical  motion  is  about  25%  greater  with 
pre-typhoon  clusters  than  with  non-developing  clusters 
when  measured  within  a large  (8°  diameter)  region  centered 
on  the  clusters.  This  difference  is  observed  primarily 
from  2-4  radius. 

2)  Pre-typhoon  cloud  clusters  are  located  in  regions  with 
mean  low-level  relative  vorticities  in  the  vicinity  of 
the  cluster  approximately  twice  as  large  as  those  observed 
with  non-developing  cloud  clusters,  i.e.  - the  ITCZ  has 
stronger  horizontal  shear. 

3)  The  cold-core  structure  in  the  500-800  mb  layer,  characteris- 
tic of  non-developing  cloud  clusters,  is  not  present  with 
pre-typhoon  cloud  clusters.  This  feature  is  further  sub- 
stantiated by  the  tangential  wind  and  relative  vorticity 


L. 


kl 


U ( m sec"') 


Fig.  19.  Zonal  wind  or  u-component  of  the  wind.  The  mean  values  of  u In 
the  r “ 0-3  region  are  plotted  for  the  pre-typhoon  (2)  and 
non-developing  (0  and  00)  clusters.  The  mean  u-component  of 
the  propagation  vector  is  also  depicted. 
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data.  The  pre-typhoon  clusters  have  a cyclonic  wind 
maximum  just  above  the  boundary  layer  as  opposed  to  a 
maximum  in  the  middle  troposphere  for  non-developing 
cloud  clusters. 

4)  The  outflow  in  the  150-250  mb  layer  is  more  anticyclonic 
with  pre-typhoon  cloud  clusters  than  with  non-developing 
clusters . 

5)  In  the  mean  non-developing  clusters  exist  in  regions 
with  larger  vertical  shear  of  the  mean  zonal  wind  al- 
though the  tropospheric  shear  of  non-developing  clusters 
may  still  be  less  tlian  5-10  m/sec.  As  a result,  the 
ventilation  integrated  through  the  troposphere  is  greater 
with  non-developing  clusters  than  with  pre-typhoon  clusters. 

This  inhibits  some  of  the  sensible  and  latent  heat  accumulations 
which  are  favorable  for  tropical  cyclone  formation. 

It  should  be  stressed  that  all  of  the  above  differences  between  pre- 
typhoon and  non-developing  cloud  clusters  were  noted  at  a time  well 
before  any  significant  intensification  of  the  pre-typhoon  cluster. 

Mean  temperature  lapse  rates  and  available  buoyant  energies  were  very 
similar.  The  observed  surface  wind  speeds  averaged  less  than  5 m sec  ^ 
in  the  immediate  area  of  both  the  non-developing  cluster  and  the  pre- 
typhoon cluster. 

Modeling  of  Genesis.  As  previously  mentioned,  all  prior  modeling 
studies  of  tropical  cyclone  genesis  have  started  with  strong,  deep 
vortices.  They  have  modeled  intensification  rather  than  genesis.  The 
transition  from  a cloud  cluster  to  a deep  closed  vortex  is  being  studied 
through  modeling  research  using  the  2 and  3-dlmenslonal  primitive  equation 
meso-scale  models  described  in  Chapter  2. 
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6.  DETERMINATION  OF  CYCLONE  INTENSITY 

It  Is  often  desirable  to  know  the  maxltnum  winds  in  a tropical  cyclone 
when  no  aircraft  reconnaissance  is  available.  At  the  present  time  the 
only  non-aircraft  data  sources  which  offer  any  promise  for  intensity 
determination  are  satellites.  One  hope  for  Improvement  in  this  forecast 
area  seems  to  be  with  establishing  a relationship  between  satellite- 
observed  eye  diameter  and  the  maximum  Intensity.  It  may  also  be  possible 
to  measure  upper  level  horizontal  temperatures  and  gradients  with  satel- 
lite or  aircraft  based  sensors  and  thereby  infer  the  low  level  wind  and 
pressure  fields.  These  topics  are  being  examined  by  Major  C.  Arnold  at 
Colorado  State  University , and  preliminary  results  are  presented  below. 

Observations  of  Eye  Diameter.  DMSP  photographs  of  57  typhoons  occur- 
ring in  the  N.W.  Pacific  from  1971-1975  were  examined.  There  were  276 
individual  observations,  and  eyes  were  visible  in  48%  of  those  photographs. 
Only  about  h of  those  eyes  were  sufficiently  well  defined  to  allow  an 
estimation  of  the  radius  of  maximum  winds  (RMW)  which  must  be  determined 
to  estimate  maximum  wind  speeds.  Therefore,  any  estimate  of  storm  in- 
tensity based  on  satellite-observed  eye  diameter  would  probably  by 
practical  about  1/4  of  the  time. 

Relationship  of  RMW  to  Maximum  Wind  Speed.  Shea  and  Gray  (1973) 
examined  the  inner-core  structure  of  hurricanes  using  NOAA  flight  data. 

They  found  a correlation  between  maximum  wind  speed  and  radius  of 
maximum  winds  (RMW)  (Fig.  21),  but  the  variability  was  so  large  that 
individual  case  predictions  would  be  tenuous.  The  results  are  the  same 
when  central  pressures  are  compared  to  satellite-observed  eye  diameters 
(Arnold,  1977).  The  central  pressure  and  maximum  winds  were  also  re- 
lated statistically  (Fig.  22),  but  variability  was  high.  It  appears 
that  the  only  way  to  estimate  storm  maximum  winds  accurately  from  eye 
dlameter/RMW  observations  is  by  using  a 2-variable  correlation  such  as 
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Fig.  21.  Variation  of  the  maximum  wind  with  RMW  for  lower  tropospheric 
data.  Tile  best  fit  curve  is  indicated  by  the  heavy  line. 
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Fig.  22.  Variation  of  the  ntaximuni  wind  with  rentral  pressure  for  all 
lower  tropospheric  data.  The  best  fit  curve  is  indicated  by 
the  heavy  line. 
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central  pressure  and  RMW  or  upper  level  (200-300  mb)  temperatures  and  I 

1 

RMW.  This  would  require  measurements  of  eye  temperatures  or  whatever  j 

parameter  besides  the  eye  diameter  was  to  be  used.  At  the  present  time  t| 

the  2-parameter  correlations  have  not  been  explored.  It  is  hoped  that 
some  method  of  relating  maximum  winds  to  satellite-observed  parameters 
can  be  developed.  Examination  of  simultaneous  satellite  and  aircraft 
data  is  required,  and  this  is  planned  for  the  immediate  future  using 
1971-1975  DMSP  and  aircraft  data  in  the  N.W.  Pacific. 

It  should  also  he  noted  that  maximum  winds  and  central  pressures 
are  not  necessarily  the  only  storm  Intensity  information  of  Interest. 

The  horizontal  extent  of  strong  winds  is  also  crucially  Important.  It 
is  unclear  how  this  parameter  could  he  determined  from  observations  of 
inner-core  characteristics. 

Relationship  of  Intensity  to  Convection.  The  current  intensity  of  a 
tropical  cyclone  is  the  cumulative  result  of  all  of  its  previous  warming. 

The  total  amount  of  deep  convection  at  any  time  Is  not  well  related  to  a 
storm's  central  pressure  or  maximum  winds.  Due  to  relatively  large  inter- 
storm  variabilities  and  diurnal  variations  in  cloudiness,  intensity 
cannot  be  predicted  well  in  terms  of  the  amount  of  deep  convection  observed. 

There  are  some  relationships  between  the  organization  of  the  convection 
and  storm  intensity  as  reported  by  Dvorak  (1975)  and  Fett  (1966).  Banding 
and  symmetry  of  convection  are  related  to  storm  intensity.  Research  into 
this  latter  type  of  intensity  determination  using  satellite  data  is 
planned  for  the  Immediate  future  at  Colorado  State  University.  It  Is 


hoped  that  by  combining  rawlnsonde  observations  with  the  DMSP  data,  more 
insight  into  the  relationships  between  convective  organization  and  inner- 
storm  strength  can  be  obtained. 


Determination  of  Intensity  from  Upper  Level  Temperatures.  Current 


research  indicates  that  surface  pressures  may  be  estimated  from  upper 
tropospheric  temperatures  (200-300  mb).  If  this  finding  is  correct, 
observations  of  temperatures  and  temperature  gradients  at  those  levels 
would  allow  determination  of  the  low  level  pressure  and  wind  fields. 

Such  temperature  measurements  could  be  made  from  direct  aircraft  obser- 
vations or  from  remote  microwave  sensors  on  satellites  or  high  flying 
aircraft.  This  is  a topic  of  current  research  at  CSU. 

Summary . It  may  be  possible  to  estimate  maximum  storm  intensity 
from  satellite  observations  of  tropical  cyclone  eye  characteristics  at 
least  1/4  of  the  time.  However,  due  to  the  extreme  variability  of 
individual  storms,  it  appears  that  observation  of  2 or  more  parameters 
is  required. 

There  are  no  significant  correlations  between  amounts  of  convection 
and  storm  intensity.  Research  into  the  relationships  between  convective 
organization  and  storm  intensity  are  planned. 

It  seems  that  storm  intensity  may  be  determined  best  from  observations 
of  upper  tropospheric  temperatures  using  satellite  or  aircraft  based 
sensors.  Research  on  this  topic  is  underway. 
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7 . CYCLONE  STRUCTURE 


The  recent  N.W.  Pacific  typhoon  rawlnsonde  composite  study  (Frank,  ; 

1976)  revealed  many  features  of  the  structure  and  energetics  of  tropical  1 

cyclones.  Some  of  the  findings  which  are  of  particular  Interest  to 
tropical  cyclone  modelers  are  discussed  below. 

Axlsymmetry . Most  current  tropical  storm  models  are  axlsymmetric 
(Kasahara,  1961;  Ogura,  1964  Yamasaki,  1968;  Ooyama,  1969;  Rosenthal, 

1970;  Sundquist,  1970;  and  Carrier  eit  ^. , 1971)  although  some  3- 
dimenslonal  models  are  operational  (Anthes  al. . 1971;  Kurlhara  and 
Tuleya,  1974;  Mathur,  1975).  Many  features  of  the  observed  storm  struc- 
ture are  relatively  symmetrical,  and  horizontal  eddy  fluxes  of  sensible 
heat  and  water  vapor  are  not  too  large.  However,  angular  momentum  and 
kinetic  energy  analyses  show  very  large  eddy  fluxes  of  these  quantities 
at  outer  radii  (Figs.  23  and  24).  Unless  these  fluxes  can  be  parameter- 
ized, It  is  doubtful  whether  an  axisymmetric  model  can  be  considered 
realistic.  Some  of  the  observed  similarities  between  the  circulations 
of  such  models  and  observed  tropical  storms  are  probably  fortuitous. 

The  eddy  generation,  dissipation,  and  export  of  kinetic  energy  are  all 
large  terms  which  tend  to  cancel,  but  this  does  not  justify  their  neglect. 

Eddy  momentum  fluxes  are  required  to  obtain  realistic  momentum  balance 
In  the  outflow  layer  at  large  radii  and  to  compensate  for  standing  eddies 
in  the  Coriolis  torque  term.  Although  many  features  of  the  tropical 
cyclone's  Inner  circulation  may  be  modeled  adequately  In  two  dimensions, 
the  authors  do  not  believe  that  such  models  can  yield  realistic  simula- 
tions of  the  outer  cyclone  regions. 
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Fig.  23.  Total  transport  of  relative  angular  momentum  (divided  by 
radius)  and  the  transport  by  the  mean  circulation  (dashed 
line)  at  r = 8 and  r = 10  . Negative  numbers  denote  import 
of  momentum. 
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Fig.  24.  Total  transport  of  kinetic  energy  (times  2)  and  the  transport 
by  the  mean  circulation  (dashed  line)  at  r - 8*^  and  r * 10*^. 
Positive  numbers  denote  export. 
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Observed  Structure.  Several  features  of  the  large  scale  structure 
of  tropical  cyclones  are  particularly  important.  Although  Inflow  at  2^ 
radius  is  virtually  all  below  800  mb,  deep  Inflow  exists  from  4°  outward 
(Fig.  25)-  This  is  not  simulated  in  most  models  which  generally  depend 
on  low  level  frictional  convergence  alone  to  force  their  radial  circula- 
tions. It  appears  that  most  of  the  middle  level  inflow  subsides  into 
the  boundary  layer  before  being  cycled  upward  in  clouds.  The  exact  cause 
of  the  mid-level  inflow  is  not  clear,  but  this  phenomenon  has  a signifi- 
cant effect  on  the  overall  storm  moist  static  energy  budget.  The  magni- 
tude of  this  Inflow  casts  doubt  on  the  use  of  CISK  as  the  exclusive 
forcing  mechanism  of  tropical  cyclones.  The  correlation  between  vorticity 
and  convergence  in  the  boundary  layer  is  poor,  and  this  also  suggests 
that  a reevaluation  of  the  CTSK  theory  is  needed. 

Soundings  (Figs.  26  a-h.  Table  6)  show  more  conditional  instability 
at  inner  radii  than  usually  is  assumed,  and  mid-level  humidities  are 
higher  than  expected.  The  inner  regions  show  substantial  potential  for 
buoyant  convection  except  in  the  eye.  The  available  data  do  not  show 
constant  vertical  profiles  of  h (or  6^)  or  moist  adiabatic  lapse  rates 
at  inner  radii.  Soundings  inside  r = 0.'^°  are  based  on  flight  data 
from  Gray  and  Shea  (1973). 

The  scale  of  the  circulation  is  very  large.  Even  at  1A”  radius 
there  is  substantial  mean  radial  flow.  Upper  level  mean  and  eddy  radial 
fluxes  of  momentum  and  kinetic  energy  are  large  at  outer  radii.  Although 
the  outer  circulation  may  be  independent  of  the  inner  core  structure, 
the  common  assumption  of  a closed  boundary  at  1000  km  radius  or  so,  does 
not  seem  justified. 
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Cunulus  Parameterization.  All  tropical  cyclone  models  must  use 
some  sort  of  cumulus  parameterization  scheme  to  simulate  the  latent 
heating  processes.  The  results  of  this  study  indicate  substantial 
convection  at  all  radii  and  large  amounts  of  vertical  mass  recycling 
everywhere  beyond  the  eyewall  region,  especially  in  the  lower  troposphere. 
Figure  27  shows  the  total  upward  motion,  mean  vortical  motion  and  re- 
quired subsidence  in  the  2-4°  outer  convective  area  based  on  a moist 
static  energy  budget  analysis.  Other  results  suggest  that  parameteri- 
zation schemes  should  include  the  effects  of  Cb  overshoot  cooling  and 
extra-cloud  subsidence  at  values  of  h higher  than  the  observed  h.  Deep 
convection  seems  well  related  to  low  level  (sfc  to  800  mb)  convergence 
as  noted  by  Cho  and  Ogura  (1974),  while  shallow  clouds  are  probably 
more  closely  related  to  low  level  winds,  turbulence  and  humidities. 

The  effects  of  ice  phase  transitions  are  probably  negligible  except  in 
areas  of  extremely  Intense  Cb  convection.  Spectral  cloud  parameteriza- 
tion schemes  may  have  to  adjust  the  types  of  clouds  used  according  to 
the  type  of  organization  of  the  deep  convection.  Vertical  momentum 
fluxes  by  clouds  should  be  parameterized  as  discussed  below. 

Sea  Surface  Fluxes.  Warm  tropical  ocean  waters  are  Important 
energy  sources  for  tropical  cyclones,  Palm^n,  1948;  1957;  Gray,  1975a. 

The  storms  characteristically  weaken  as  they  move  over  colder  v’aters 
(Brand,  1971). 

The  sea  surface  to  air  fluxes  of  moist  static  energy  (h)  found 
in  the  recent  study  of  Frank  (1976)  are  smaller  than  those  reported  in 
most  previous  empirical  studies  (Table  7).  The  results  of  the  current 
study  should  be  the  best  estimate  due  to  the  very  large  size  of  the  data 
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extra-cloud  subsidence  (oi) . The  effect  of  Including  ice  phase 
transitions  is  shown. 


TABLE  7 

2 

Sea  Surface  to  Air  Fluxes  of  Moist  Static  Energy  (h)  (cal/cm  day) 


Study 

h Fluxes 

Region 

Present  Study 

1470 

(0-80  km) 

820 

(80-220  km) 

Daisy  (Rlehl  & Malkus) 

August  25 

1019 

(0-130  km) 

August  27 

1675 

(0-130  km) 

Hilda  (Hawkins  & Rubsam) 

-^2560 

(16-130  km) 

Hilda  (Lelpper) 

'V.4150 

('\0-240  km) 

Inez  (Hawkins  & Imbembo) 

'v-4018 

(18-90  km) 

(Malkus  & Rlehl) 

3140 

(30-90  km) 
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A constant  sensible  heat/evaporation  coefficient  of  1.5  x 10 
applied  to  the  900  mb  winds  seems  adequate  to  compute  fluxes  of  these 
quantities  from  classical  bulk  aerodynamic  formulas.  A variable  drag 
coefficient  seems  to  yield  better  estimates  of  mt-nnentum  and  kinetic 
energy  fluxes  at  the  sea  surface. 

Radiation . Diurnal  variations  in  tropical  storm  convection  and 
temperatures  are  observed  indicating  that  cloud/cloud  free  radiation 
differences  should  be  incorporated  into  tropical  cyclone  models.  Pre- 
cipitation is  maximum  in  the  late  morning  and  minimum  in  the  evening,  but 
the  variation  is  only  about  15%.  Middle-level  temperatures  are  warmest 
during  the  day.  The  radlatlonal  forcing  occurs  day  and  night  and  may 
well  be  a factor  in  such  Important  storm  features  as  the  middle  level 
inflow.  Since  the  dominant  factor  in  the  radiation  budget  is  the 
existence  of  dense  cirrus  cover,  the  extent  of  the  cirrus  shield  must 
be  predicted.  The  Importance  of  radlatlonal  forcing  is  also  stressed 
in  observational  studies  by  Jacobson  and  Cray  (1976)  and  Foltz  (1976). 

Angular  Momentum.  In  addition  to  the  Importance  of  horizontal 
eddy  transports  of  momentum  mentioned  above,  it  must  be  noted  that 
the  Coriolis  torque  term  (fV^r)  does  not  integrate  to  zero  in  a region 
containing  persistent  northerly  or  southerly  flow.  At  large  radii  this 
may  become  a major  term  in  the  angular  momemtum  budget. 

Vertical  eddy  fluxes  of  momentum  by  cumulus  clouds  are  required 
to  explain  the  observed  circulation.  They  should  be  included  in  cumulus 
parameterization  schemes  wherever  vertical  wind  shears  are  significant. 

As  a first  approximation,  angular  momentum  can  be  treated  as  a conserved 
quantity  in  clouds. 

Kinetic  Energy.  The  use  of  mean  circulation  quantities  alone 
gives  a poor  facsimile  of  the  true  kinetic  energy  budget.  The  sea  surface 
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dissipation,  mean  circulation  transport,  and  mean  circulation  generation 
terms  tend  to  balance.  This  has  often  led  to  the  neglect  of  eddy 
processes.  However,  eddy  fluxes  are  large,  and  Internal  dissipation  and 
eddy  generation  are  substantial  - perhaps  dominant.  The  large  scale 
storm  circulation  Is  a strong  source  of  kinetic  energy  exporting  about 
2 

A. 8 watts/m  averaged  over  the  0-10  radius  region.  This  Indicates  that 
tropical  cyclones  may  play  a significant  role  in  the  tropical  general 
circulation  during  certain  seasons. 

Water  Budget . Due  to  the  basic  pattern  of  lower  and  middle  tropo- 
spheric inflow  and  upper  tropospheric  outflow,  tropical  cyclones  converge 
large  quantities  of  water  vapor.  This  phenomenon  is  enhanced  by  above 
average  surface  evaporation  of  sea  water  due  to  high  surface  winds.  A 
small  part  of  this  water  vapor  convergence  goes  into  Increased  humidity 
levels  in  the  inner  storm,  but  most  falls  out  as  precipitation. 

The  vertically  integrated  water  budget  of  the  composite  steady  state 
t5rphoon  is  shown  in  Fig.  28.  Horizontal  convergence  of  q was  observed, 
and  precipitation  was  estimated  from  rainfall  studies  by  Miller  (1958) 
and  Gray  et  al . (1975).  Evaporation  was  computed  as  a residual.  The 
results  seem  realistic  except  for  an  apparent  overestimate  of  evaporation 
from  4-6°  and  an  underestimate  from  6-8°  - probably  the  result  of  a slight 
error  in  the  r = 6°  radial  wind  profile. 

It  is  interesting  to  note  that  while  about  75%  of  the  observed 
0-2°  precipitation  results  from  horizontal  convergence  of  water  vapor, 
nearly  65%  of  the  0-6°  precipitation  may  be  attributed  to  evaporation 
within  that  area.  The  0-6°  evaporation  is  about  double  mean  tropical 
values,  primarily  due  to  high  winds.  Hence,  nearly  1/3  of  the  0-6° 
precipitation  results  from  this  anomalous  evaporation.  Tropical  cyclones 
extract  large  mounts  of  extra  energy  from  the  sea. 
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I 

WATER  BUDGET  (cm/day) 

CONVERGENCE  ,2  -q?  q.I  -O.I  -O.I 

OF  q 

PRECIPITATION  9.0  2.3  0.7  0.7  0.7  0.7 

REQUIRED 

EVAPORATION  I 1 2.3  I t l.l  I | 1.4  I j 0.6  I 1 0.8  I |0.8 

(T  2°^  4®^  6^^  8°^  10^^  12° 

Fig.  28.  Water  budget  for  0-12°  region.  Evaporation  values  shown 

are  the  amounts  required  to  achieve  balance.  j 

i 
i 

Typhoon  Rainfall.  N.W.  Pacific  rainfall  data  from  small  Island 
stations  were  analyzed  to  determine  the  mean  rainfall  rates  by  quadrant 
for  the  0-4°  radius  area  (Fig.  29).  No  corrections  were  made  for 
orographic  effects  or  measurement  errors.  The  values  shown  are  for  the 
0-2°  band  and  the  2-4°  band  for  the  four  quadrants  (the  grid  is  aligned 
with  the  direction  of  storm  motion  toward  the  top  of  the  page) . Also 
shown  are  the  radial  band  averages  for  all  four  quadrants.  No  striking 
asymmetries  are  observed,  but  rainfall  does  seem  to  be  highest  In  the  right 
rear  quadrant.  This  finding  contradicts  earlier  estimates  of  strong  rain- 
fall maxima  In  the  right  front  quadrant  (Cline,  1962;  Schoner,  1957). 

Since  most  previous  tropical  cyclone  rainfall  studies  have  been  performed 
for  storms  at  or  near  continental  landfall.  It  Is  possible  that  the 
previously  observed  right  front  quadrant  maxima  are  orographic  In  origin 
due  to  onshore  flow.  Under  such  conditions  supressed  convection  would 
be  expected  In  the  front  left  quadrant  where  off-shore  flow  prevails. 


6.7 

1.2 

-0.7 

0.1 

-0.1 

-0.1 

9.0 

2.3 

0.7 

0.7 

0.7 

0.7 

1 2.3 

.4^ 

j 0.6 
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COMBINED  HOURLY  B JAPANESE  DATA 


(f-<Bon«  4987nn  IIOS  cm  S.30 

Fig.  29.  Mean  precipitation  around  tropical  cyclone  based  on  hourly 
data  from  9 Island  stations  and  12  hourly-data  from  4 
Japanese  Island  stations. 


Dunn  and  Miller  (1964),  have  proposed  that  the  precipitation  pattern  of 
low  latitude  storms  Is  symmetrical  with  a right  front  quadrant  maximum 
developing  as  the  storms  move  poleward. 

The  hourly  rainfall  data  were  also  used  to  determine  the  frequency 
and  Intensity  of  rainfall  occurrences  during  typhoon  passage  near  an 
Island  station.  Table  8 shows  the  percentages  of  the  total  time  In 
which  various  rainfall  rates  were  recorded  for  each  radial  band. 

It  Is  Interesting  to  note  the  relatively  small  percentage  of  the 
time  that  moderate  to  heavy  rains  .1  In/hour)  fell.  Such  rainfall 
rates  were  recorded  only  35%  of  the  time  In  the  highly  active  0-2° 
region  and  but  17%  of  the  time  In  the  2-4°  band.  Outside  the  4°  radius 
rainfalls  of  this  intensity  were  not  often  recorded  (3-6%  of  the  time). 
Substantial  amounts  of  light  rain  (.01-. 10  in/hour)  were  recorded  Inside 
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TABLE  8 

Rainfall  Rates 


Rainfall  Rate 

Percent 

of  Tot 

al  Storm 

Passage 

Time 

in /hour 

0-2° 

2-4° 

^6° 

1-  ! 
too 

0 

^-10° 

10-12° 

12-14° 

0-trace 

30 

59 

80 

87 

89 

89 

90 

.01-. 1 

35 

24 

14 

10 

7 

7 

7 

.1-.3 

19 

11 

4 

2 

3 

2 

2 

.3-. 9 

13 

5 

2 

1 

1 

2 

1 

>.9 

3 

1 

<1 

<1 

<1 

<1 

<1 

6 . It  is  possible  that  some  of  this  rain  was  due  to  orographic  effects 
resulting  from  elevated  terrain,  strong  winds,  and  very  high  humidities. 
The  above  average  rainfall  recorded  in  the  4-6*^  band  seems  to  result  pri- 
marily from  an  above  normal  incidence  of  light  rainfall  and  slightly  above 
average  moderate  rainfall  (.1-.3  In/hour).  Since  this  area  is  observed 
to  be  relatively  cloud  free,  except  for  occasional  outer  ralnbands,  at 
least  part  of  the  observed  positive  anomaly  in  precipitation  may  have 
been  orographic  in  origin. 

The  results  of  this  rainfall  study  Indicate  that  typhoon  precipi- 
tation is  primarily  concentrated  in  localized  heavy  convective  regions. 

The  incidence  of  2-4°  band  rainfall  is  only  slightly  greater  than  was 
found  for  N.W.  Pacific  cloud  clusters  by  Ruprecht  and  Gray  (1976).  The 
inner  core  (0-2°)  area  shows  a much  larger  incidence  of  rainy  episodes 
than  is  found  in  cloud  clusters  as  expected.  The  frequencies  of  moderate 
to  heavy  rainfall  (>.l  In/hr)  agree  relatively  well  with  the  convective 
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cloud  covers  inside  200  miles  of  *^5%  Cb  and  '^2SX  total  cumulus  cloud  cover 


reported  by  Malkus  et  al . (1961)  and  Gentry  (19b4). 
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8 . Sl'MMARY 

A number  of  new  findings  regarding  the  structure,  dynamics  and 
behavior  of  tropical  cyclones  have  been  summarized.  However,  much  addi- 
tional observational  and  modeling  work  remains.  An  enormous  data  set 
using  many  types  of  observations  from  various  regions  has  been  assembled 
at  CSU . More  data  is  continually  being  acquired.  This  is  the  largest 
tropical  cyclone  data  sample  assembled.  The  results  of  this  paper  are 
in  many  ways  only  preliminary  findings.  The  researchers  are  confident 
that  current  and  future  research  efforts  will  significantly  Improve 
tropical  cyclone  forecasting  skills  and  Improve  our  overall  understand- 
ing of  tropical  weather  systems. 
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